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Abstract 

The problem of the orientational ordering transition for lattice-gas mod- 
els of liquid crystals is discussed in the low- dimensional case d = 1,2. For 
isotropic short-range interactions, orientational long-range order at finite 
temperature is excluded for any packing of molecules on the lattice Z d ; on the 
other hand, for reflection-positive long-range isotropic interactions, we prove 
existence of an orientational ordering transition for high packing (// > /^o) 
and low temperatures (/3 > /3 C (//)). 
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The problem of long-range orientational order in models mimicking liquid 
crystalline (usually nematic) behaviour has been discussed at the rigorous 
level and in the framework of different interaction models, see, e.g, Ref. 

11,11- 

In the following, we shall be considering cylindrically symmetric molecules, 

whose centres of mass belong to a d— dimensional space (x G Z d , lattice 

models, or x G R d , continuum models), and whose orientations are defined 

by m— component unit vectors u G S^" 1-1 ) C R m , m > 2. 

For low-dimensional lattice models (i.e. dimension d — 1,2) with isotropic 

short-range interactions, the Mermin- Wagner theorem entails absence of an 

orientational ordering transition taking place at finite temperature; a similar 

result was proven by Romerio for continuum fluids |l|. On the other hand, 

a two-dimensional lattice model with anisotropic interactions restricted to 

nearest neighbours can produce an ordering transition (also of the nematic 

type), see, e.g. Ref. §, Ref. @, Ch. 9, Example 9.22(2). 

The proof of the existence of orientational order for continuum liquids is a 

rather complicated problem. In this connection we would like to mention 

Ref. J7|, where the existence of a ferromagnetic phase transition is proven 

for a continuum fluid (d > 2) of classical particles carrying Ising-like spins 

and having suitable magnetic and non-magnetic interactions. The proof is 

based on the FKG- and GHS-inequalities, which make it possible to find a 

lower bound for the magnetization of the fluid in terms of that of a suitable 

spin system on a lattice; no such inequalities are known for m > 3. 

The aim of the present note is to prove the existence of orientational ordering 

(at finite temperature) in low- dimensional lattice-gas models (d = 1, 2) M 

with isotropic long-range interactions of nematic symmetry, which therefore 

do not satisfy the hypotheses of the Romerio theorem M. 



In order to make contact with (nematic) liquid crystals, and following the 
standard phenomenology ||, we consider centrosymmetric molecules (sym- 



metry -Doo/i), and define the matrix 

Q^ = («V--y, a, (3 = 1,2,... m; (1) 

m 

then translationally invariant interactions between molecules can be ex- 
pressed as 

V xy = -a(\x-y\)Tr(Q x -Q y ) + g{\x-y\), x,y G R d . (2) 

The Hamiltonian of the lattice-gas caricature of the nematic liquid crystal 
has the form Q: 

H A (h, Q) = -- Y^ n x n y J xy Tr(Q x ■ Q y ) + - ^ n x n y I xy - \i ^ n s , (3) 

(»>») (a:,!/) x e A 

where (x, ?/) = {x,y 6 A:i/ y}. In this context we can assume a(0) = 
(/(O) = 0. 

So, molecules live on the sites of the cubic sublattice k d Z d with periodic 
boundary conditions, i.e. 

J xy = J2 a(\x-z\), I xy = J2 9{\x-z\). 

{zeZ d :z=y( mod A)} {z£Z d :z=y( mod A)} 

A configuration is specified by a set of occupation numbers {n x — 0, 1; x G 
Z d } and, for d\\ {x E Z d : n x = 1} , by the configuration Q x of the molecule 
at x G Z d . The corresponding finite-volume Gibbs state is defined by 

(f} A ((3, f ,) = [E A ((3^)}- 1 
E / II (dv(Q x )T*exp[-[3H A (n, Q)]f{n, Q), (4) 

{n x =0,l; xGA} xeA 

where Sa is the partition function, dv is the 0(m)— invariant probability 
measure induced by the Haar measure on the unit sphere in R m , see Eq. ([[]), 
and n = {n x ; x G A}, Q = {Q x ; x G A}. 

The chemical potential \x governs the concentration (mean density) p(/3, /j) = 
(n x ) A , of molecules on the lattice at the temperature f3~ l . The interaction 
term J xy Tr(Q x ■ Q y ) involves both positional and orientational degrees of 



freedom, and possibly produces orientational order, whereas I xy mimics a 
direct interaction between them (positional order). 

Let a(\x — y\) denote a short-range interaction (SR, i.e. finite-ranged or 
decreasing at least exponentially), or a long-range one with asymptotic be- 
haviour a(\x — y\) ~ \x — y\~( d+a \ for \x — y\ — > oo, a > 0. Then, in the limit 
q — > 0, the lattice Fourier transform of J xy is given by 

J(0) - J(?) * { c\q\\ a>2,orSR (5) 

where 

q e A * = { g « = ^m,, m a = 0,±1,...,±(^ - 1),+^, a = 
1,2,. ..rf}. 

Then, since this interaction is 0(m)— invariant, we have the following Propo- 
sition a la Mermin- Wagner, due to ||. 
Proposition 2.1 

Let d — 1, 2, a > 2 (or let a(|x— y|) define a 57? interaction), and let #(|:e— J/|) 
be such that 

Yl \9(\x-y\)\ \x-y\ 2 < +00. 

y£Z d 

Then there is no orientational order at finite temperature, i.e. 

P(P^)= lim lim (Tr{D x ■ D y )) A = 0, (6) 

|x— 1/|— *oo AyZ d 

Here D x = n x ■ Q x . 
Proof 

For any fixed configuration n the limiting Gibbs state (—)(/3,n) is 
0(m)— invariant. In this case, the proof developed in || carries through 
verbatim for the 0(m)— invariant interaction defined by (0) and correspond- 
ing to U xy (Q x , Q y ) in Pfister's notation 0. Since the conditions both on I xy 
and J xy guarantee for the Hamiltonian (|3p the equivalence of ensembles (see 



Ref. [II], [TT|]), one gets that the grand-canonical Gibbs state 

<J)(P,H)= I dK(n,n)(f)(P,h) (7) 



is also 0(m)— invariant. Here dK(fi,n) is the Kac-transformation kernel 



K], [llj. For any pure state (—)(/?, /i), we get 



lim (Tr(D x ■ D y )) = (TrD x ) 2 = 0, (8) 

where the last equality follows from the 0(m)— invariance of the state (|7|). □ 
The presence of "holes" (n x = 0) evidently disfavours the orientational order 
on lim A | Z d M(n). It is natural to guess that the order parameter P(/3,n), 
and hence P(j3, /i) are bounded by the order parameter of the lattice without 
"defects", i.e. {n x — 1; x & Z d }, which corresponds to /i — > +oo. One can 
easily check this for particular cases of the lattice-gas Ising system (where 
Tr(Q x ■ Q y ) is substituted by t x ■ r y , r x ^ y = ±1), or for the lattice-gas of 
plane rotators, due to the GHS inequalities. Hence, for these systems, the 
result (f|) is a simple consequence of "pure system domination" : P(/3, fi) < 
P(/3,fi = +oo), and of the Mermin- Wagner theorem for the regular lattice 
Z d , d=l,2. 

This result can be futher strengthened: for d = 1 and 1 < a < 2, the absence 
of orientational order has also been proven for the lattice model(s) [Ref. |J, 



Ch. 9, Comment 9.34; Theorem 14' in Ref. Il2|; Refs. I3L HI . Notice that 



Romerio's result fllfl holds for a > 2 but for a continuum model. 



In contrast to the "no-go" Proposition 2.1, the proof of existence of orienta- 
tional order at finite temperature (P(/9,/x) > 0) is a more delicate task. We 
are able to do this for long-range interactions which are reflection-positive 
with respect to reflections in planes without sites |12|, |i~5f . 
From now on, let g{\x—y\) = and let a(\x— y\) correspond to the long-range 
interaction 

a(\x - y\) = b\x - y\- {d+a) , b > (9) 

with < o < d. First we note that, for (|9|) the asymptotic form of the 
Fourier transform (|5]) excludes a a "no-go" theorem a la Mermin- Wagner. 
Proposition 3.1 

For d = 1, 2 there exists a fio and, for every \i > /i there exists a (3 C {^) such 
that P((3,fi) > for j3 > /3 c (n). This means that every limiting Gibbs state 
(— ) = lim A | Z d(— )a has long-range orientational order. 



Proof 

The Proof is an adaptation to our case (|l|) of the line of reasoning developed in 

[|j] for the case of the general matrix order parameter. By (|I|) one has TrQ x = 

0, i.e. Tr{D x ) A = (TrD x ) A = 0. Then , by the 0(m)— invariance of the state 

(— )a, one can deduce that (D x )^ = 0. Hence, P(j3,fi) > Tr((D x ) 2 ) = will 

really mean existence of long-range orientational order. 

Using the translational invariance of (— }a> we get 

c A = \A\- 2 J2(Tr(D x D y )) A = Tr(D 2 x ) A -\A\- 1 ]T Tr(D p ■ LL p ) A , (10) 

(x,y) peA*\{0} 

where 

D p = |A|- 1 / 2 ExeAexp(-zpa;) J D,. 

Thus, in order to prove that P([3,fi) > 0, it suffices to verify that 

c(3,u) = lim inf c A > 0. (11) 

K\Z d 

The upper bound on Tr(D p D_ p ) A for p e A* \ {0} results from the Infrared 
Bound P, |15| 



~ — con if 

(Tr(D p D_ p )) A < ~ \ P \-, forp^O (12) 

/3[J{0) - J{p)\ 

which holds true for the reflection-positive interaction @ , due to the chess- 
board estimate proving the gaussian domination; see Refs. [^, [LlJ for a 
review. Therefore, the sum over A* \ {0} divided by |A| (see fllCf)) can be 
estimated in the limit A ] Z d from above by the integral 

h,M = ^ / d d p[m - aipT 1 <oc,0<a<d. (13) 

In order to get the lower bound of the first term of (0), Tr(D 2 ) A = 
(nl)\Tr(Ql) = I2 — I^^a, we again use the chessboard estimate |~2], [Lj| : 



)A<[(ri( i -^))A] 1/|A - (14) 



n ~, 

yeA 



Then one gets 



^A^i-Kna-^))^ 1 ^- (15) 



yeA 



According to (|J) the last term in ([15]) is equal to [Za(/3,/x)] _1 /' a L To get a 
lower bound on the partition function Z\(fi, //), we choose a D e swpp 1/ and 
a neighbourhood jV e of D such that 

Tr(D l -D 2 ) > (l-e)Tr(D 2 ) > 0, (16) 

for some e > and for every Di, D 2 & J\f e . Then for configurations 

A^ = {D x -.D x eK,xeA}, 

we get by (|16"D that {n x = 1; i6A} and 

- #a(ti, Q) > (1 - e)TrQ 2 £ J^ + M |A| = [(1 - e)TrQ 2 a(0) + M ] ■ | A| (17) 

Therefore, upon restricting the integration in the partition function to con- 
figurations A^ e \ we obtain 

\Z A {[5^)}- 1 ^ < HK)T l exp{-/3[(l -e)TrQ 2 a(0) + ft]}, (18) 

and, as a consequence of (|15|) , one gets 



Tr(Dl) A > jj[l - ^yexp{-/5[(l - e)TrQ 2 a(0) + //]}] = L(/?,/i) (19) 

Combining (p~3| ) with ([19]), we obtain the following estimate for (|TT]): 

c%n)>L%n)-I^). (20) 

Now one immediately sees that for 
/i>/i = -(l-e)Tr(Q 2 )a(0), 
there exists a /3 c (/x) such that 
L(p c { t i),n) = I d , a {p), 

where < cr < 1, d—1 and < cr < 2, d = 2. Hence, for /3 
beta c (fj, > /io) we get c(/5,/i) > 0. According to Eqs. QIDD, (Q) 5 this means 
that P(/3, /i) > in the named (/3, /i) domain, i.e. in any pure limiting Gibbs 
state (D x ) = (Q x ) 7^ 0: thus the 0(m) symmetry is broken, which means 
long-range orientational order. □ 

Notice that the sign of a(\x — y\) plays an important role in the existence 
theorem, but not in its absence counterpart. 

Moreover, the following corollaries can be obtained, by the same line of rea- 
soning as in the preceding propositions. 

6 



Corollary 1. 

Consider the ferromagnetic (lattice gas) counterpart of the present model, 

i.e. whose orientation-dependent two-body interaction reads 

V xy = -a{\x-y\){u x -u y ). (21) 

Then, under the same hypotheses on a(\x — y\) and g(\x — y\, one can prove 

absence or existence of an ordering transition, respectively. 

Corollary 2. 

Let m = 3, and let the orientation-dependent two-body interaction have the 

general form 

V xy = -a(\x-y\)P L (u x -u y ), (22) 

for an arbitrary positive integer L; owing to the addition theorem for spherical 
harmonics, the Legendre polynomial can be written as 

P L (u x ■ u y ) = u L Tr(\N x ^ • WW), (23) 

where W' ' denote real multipole tensors of rank L constructed in terms of 
components of u x and u y , respectively, and ujl is an appropriate positive 
normalization factor, so that Tr(\N x ' ■ \N X ') = 1/ujl- Then, under the same 
hypotheses on a(\x — y\) and g(\x — y\), one can prove absence or existence of 
an ordering transition, respectively, for arbitrary L. Even values of L define 
lattice-gas models of nematic liquid crystals, and the result proven in Ref. 
[|j] and for d = 3 can be similarly generalized. 

In conclusion, we would like to point out that the present note yields a partial 
answer to questions proposed in [1[] as open problems. We have shown that 
long-range interactions violating conditions (2.5) of Ref. |IJ (cf. @ ) can 
produce long-range orientational order in low-dimensional lattice models of 
nematic liquid crystals. The problem of the possible existence of long-range 
orientational order for continuum models remains open. 
Another open problem which can be formulated in the framework of this 
lattice model (|3|) concerns existence of long-range positional order, i.e. 
lim| x _j / |_ >+00 [(n :!; • n y ) - (n x ){n y )] > 0, 

when g{\x — y\) > 0, repulsion which crudely mimics excluded- volume effects 
for nematic molecules. 



Acknowledgments 

The present text reached its final form during S.R.'s visit to CPT-Luminy; 
the French CNRS is acknowledged for hospitality and financial support. 



References 

[1] M. Romerio, J. Math. Phys. 19 (1978) 802. 

[2] O. J. Heilmann and E. H. Lieb, 

Commun. Math. Phys. 25 (1972) 190; 
J. Stat. Phys. 20 (1979) 679. 

[3] D. B. Abraham and O. J. Heilmann, J. Phys. A 13 (1980) 1051. 

[4] N. Angelescu and V. A. Zagrebnov, J. Phys. A 15 (1982) L639. 

[5] S. B. Shlosman, Russian Math. Surveys 41 (1986) 83. 

[6] H. O. Georgii, 

Gibbs Measures and Phase Transitions 
(De Gruyter, Berlin, 1988). 

[7] Ch. Gruber and R. B. Griffiths, Physica A 138 (1986) 220. 

[8] P. G. de Gennes and J. Prost, 

The Physics of Liquid Crystals, second edition 
(Oxford University Press, Oxford, 1993). 

[9] C.-E. Pfistcr, Commun. Math. Phys. 79 (1981) 181. 

[10] J. Lewis, C. E. Pfister and W. G. Sullivan, 

Large Deviations and Thermodynamic Formalism: a New Proof 
of the Equivalence of Ensembles, in 

On Three Levels: Micro-, Meso- and Macro- Approaches in Physics, 
eds. M. Fannes, C. Maes and A. Verbeure, 

NATO ASI Series B, Plenum Press (London-New York, 1994); 

proceedings of the Conference held in Leuven (July 19-23, 1993). 



[11] J. Lewis, C.-E. Pfister and W. G. Sullivan, 

The Equivalence of Ensembles for Lattice Systems: Some 
Examples and a Counterexample, 

J. Stat. Phys., in press. 

[12] J. Frdhlich, Bull. Am. Math. Soc. 84 (1978) 165. 

[13] G. Nenciu, Phys. Lett. A 34 (1971) 422. 

[14] P. A. Vuillermot and M. V. Romerio, Commun. Math. Phys. 41 (1975) 
281. 

[15] J. Frdhlich, R. Israel, E. H. Lieb, and B. Simon, Commun. Math. Phys. 
62 (1978) 1. 

[16] Ya. G. Sinai, 

Theory of Phase Transitions: Rigorous Results 
(Pergamon Press, Oxford- New York, 1982). 



10 



